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Abstract 

As a consequence of dynamical chiral symmetry breaking the pion Bethe- 
Salpeter amplitude necessarily contains terms proportional to 75 7 • P and 
75 7 • k, where k is the relative and P the total momentum of the constituents. 
These terms are essential for the preservation of low energy theorems, such as 
the Gell-Mann-Oakes-Renner relation and those describing anomalous decays 
of the pion, and to obtaining an electromagnetic pion form factor that falls as 
1/q 2 for large q 2 , up to calculable In ^-corrections. In a simple model, which 
correlates low- and high-energy pion observables, we find q 2 F 7r (q 2 ) ~ 0.12 - 
0.19 GeV 2 for q 2 > lOGeV 2 . 

Pacs Numbers: 13.40.Gp, 14.40.Aq, 12.38.Lg, 24.85.+p 
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I. PION AS A BOUND STATE 



Understanding the pion is a key problem in strong interaction physics. As the lowest 
mass excitation in the strong interaction spectrum it must provide the long-range attraction 
in N-N potentials 0. In QCD, it is a quark-antiquark bound state whose low- and high- 
energy properties should be understandable in terms of its internal structure, and it is 
also that nearly-massless, collective excitation which is the realisation of the Goldstone 
mode associated with dynamical chiral symmetry breaking (DCSB). An explanation of these 
properties requires a melding of the study of the many body aspects of the QCD vacuum 
with the analysis of two body bound states. The Dyson- Schwinger equations (DSEs) provide 
a single, Poincare invariant framework that is well suited to this problem. 

The DSEs are a system of coupled integral equations and truncations are employed to 
define a tractable problem. In truncating the system it is straightforward to preserve the 
global symmetries of a gauge field theory and, although preserving the local symmetry 
is more difficult, progress is being made 0. The approach has been applied extensively (4j 
to the study of confinement, and to DCSB where the similarity between the ground state of 
QCD and that of a superconductor can be exploited, with the QCD gap equation realised 
as the quark DSE. It has also been employed in studying meson-meson and meson-photon 
interactions ||, heavy meson decays 0, QCD at finite temperature and density [0, and 
strong interaction contributions to weak interaction phenomena |HJ. 

Studying the pion as a bound state requires an understanding of its (fully-amputated) 
Bethe-Salpeter amplitude, which has the general form 



ri(k;P)=r j l5 



iE v (k;P)+~fPF v (k;P) (1) 



+ ~f-kk-P G^k- P) + a k^P v H^k; P) 

where {r jf }j =1 3 are the Pauli matrices. r£ satisfies the renormalised, homogeneous Bethe- 
Salpeter equation 

~K(k;P)] tu = j\xi(q;P)} sr K™( q ,k;P), (2) 

where k is the relative and P the total momentum of the quark-antiquark pair, xHq'i P) := 
S(q + )r{(q; P)S(q_), r,. . . ,u represent colour, flavour and Dirac indices, q± = q ± P/2, and 
J g A := J A (i 4 g/(27r) 4 represents mnemonically a translationally-invariant regularisation of the 
integral, with A the regularisation mass-scale: A — > oo is the last step in any calculation. In 
Eq. (H), K is the fully-amputated, renormalised quark-antiquark scattering kernel and 5* is 
the renormalised dressed-quark propagator, which is the solution of 

Sip)' 1 = Z 2 {ij -p + m hm ) + Z 1 g 2 D fM/ (p- qh»S(q)T v (q,p) , (3) 

Jq 

where D^ v {k) is the renormalised dressed-gluon propagator, r M (g;p) is the renormalised 
dressed- quark-gluon vertex and mb m (A) is the Lagrangian current-quark bare mass. In 
Eq. (H), Z\ and Zi are the renormalisation constants for the quark-gluon vertex and quark 
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wave function, and the chiral limit is obtained with mt, m (A) = 0. The solution of Eq. ([3D 
has the general form 



S(p) = -17 ■po v (p 2 ) + a s (p 2 



«7 • p A(p 2 ) + B(p 2 



(4) 



Also important in the study of the pion is the chiral-limit, axial- vector Ward-Takahashi 
identity 



- iP»T^{q- P) = T5 y + 75 y S-\q-) 



(5) 



This identity relates the renormalised, dressed-quark propagator to the renormalised axial- 
vector vertex, which satisfies 



tlx 



757^ - 



- tu 



[ X Uq;P)] sr K r t :( q ,k;P) 



(6) 



where xiJQl P) '■= S(q+)T J 5fl (q; P)S(q-). In the chiral limit the general solution of Eq. 
is 



10,11 



!%,(*; P) = y75 



l^Fn{k\ P) + 7 • kk^G R {k; P) - <v k » H R {k; P) 

p 

o r n 



(7) 



+ rl fM (k;P)+ fijgr^P) 



where: F R , G R , H R and are regular as P 2 -> 0; P^i^k-P) ~ 0(P 2 ); r w (ife;P) is 
the amplitude in Eq. (|Ij); and the residue of the pole in the axial- vector vertex is /°, the 
chiral-limit leptonic decay constant, which is obtained from: 



S iS f«Pu 



tr 



T 



J 5 l»S(q + )ri(q;P)S(q- 



(8) 



[This expression is valid for arbitrary values of the quark mass.] Now, independent of 
assumptions about the the form of K, it follows [TO] from Eqs. (51), (Bp and (|7]) that 



f°E„{k 1 0) = B (k 2 ), 
F fl (fc;0) + 2/X(A;;0) = Ao(A; 2 ), 
G R (k;0) + 2f°G„(k;0) = 2A' (k 2 ), 
H R (k;0) + 2f°H 7T {k;0) = 0, 

where A and B are the chiral limit solutions of Eq. (j^). A necessary consequence of Eqs. @- 
(1^) is that the pseudovector components F n and G n , and the pseudotensor component H n , 
are nonzero in Eq. ([!]). This corrects a misapprehension [OJ that only i?^ ^ and has 
important phenomenological consequences. 



(9) 
(10) 

(11) 
(12) 
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A. Normalisation of the pion field 



To highlight one such consequence we note that Eq. fl2|), the homogeneous Bethe-Salpeter 
equation, does not determine the normalisation of the Bethe-Salpeter amplitude. The canon- 
ical normalisation is fixed by requiring that the pion pole in the quark-antiquark scattering 
amplitude: M := K/[l — (SS)K], have unit residue. As an alternative, one can normalise 
the solution of Eq. (|2]) by requiring that E(0; 0) = B(0) in the chiral limit. In terms of the 
amplitude Q n (k; P) defined in this way, the canonical normalisation condition is 



^ 3 Nl P, = jf A jtr QUq; - p )^^Gi(Q; P)S(q 

dS{q„ ' 



(13) 



+tr 



gi(q;-P)S(q + )gi(q;P) 



dP„ 



+ 



A M 



[Xl(q; -P)] /g7 1 J [Xk(k; P)] ut , 



where Q^q; — P)* := C^ 1 Q n (—q; —P) C with C = 7274, the charge conjugation matrix, 
and X f denoting the matrix transpose of X. Equation (|U|) defines the pion normalisation 
constant, N^, which has mass-dimension one. Physical observables are expressed in terms 
ofT 7r (k;P):=(l/N w )g 7T (k;P). 

In the chiral limit, when all the amplitudes in Eq. ([!]) are retained, one obtains flCfl 



f° = N°. 

J TT TT 



(14) 



This result verifies a core assumption in chiral perturbation theory; i.e., that the pion field 
is normalised by /°, which is implicit in the expression of the chiral field as 



U(x) 



jT-n(x)/f° 



(15) 



However, Eq. ( |14|) is violated in bound state treatments of the pion that neglect the pseu- 
dovector components of the Bethe-Salpeter amplitude |TT[ .[] 



II. ANOMALOUS NEUTRAL PION DECAY 



7T 



The pseudovector components of the pion also play a special role in the anomalous 
— > 77 decay. Consider the renormalised, impulse approximation to the axial-vector- 



photon-photon (AW) amplitude:^ 

Tpnu{ki, k<i) := T p ^ u {kx, k 2 ) + T pu ^[k2i k\, 



T P p U (ki, k 2 ) 



tl£)F 



S (qi) T 5 P (q; -P) S(q 2 ) in(ga, qi2)S(q 12 ) iTZ(qu, qi) 



(16) 
(17) 



1 N 7T in Eq. (jTJ|) provides the best numerical approximation to the pion's leptonic decay constant 
in analyses that neglect the pseudovector components and employ a P-independent form for K. 



2 In our Euclidean metric: 7^ = 7^, {7^,7^} = 2 6^, and a spacelike vector, kp, has k 2 > 0. 
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where ki, k 2 are the photon momenta [k 2 = = k 2 , 2k\ ■ k2 = P 2 }, q is the loop- momentum, 
and q x :=q- k x , q 2 := q + k 2 , q := f (<?i + q 2 ), ?ia ■= q - h + k 2 . 

Here E^px,^) is the renormalised, dressed-quark-photon vertex, and it is because this 
vertex satisfies the vector Ward-Takahashi identity: 



(Pi -P2) t iiFl(p 1 ,p 2 ) = S \pi)-S 1 (p 2 ), 



(18) 



that no renormalisation constants appear explicitly in Eq. fll7|) . VUp^p^) has been much 
studied || and, although its exact form remains unknown, its robust qualitative features 
have been elucidated so that a phenomenologically efficacious Ansatz has emerged |L3 



iTl(p, q) := z£ A (p 2 , q 2 ) 7m + (p + g) M [|i 7 • (p + g) A A (p 2 , q 2 ) + A B (p 2 , g 2 ; 
*ttfrf)-=\W) + ftf)\, 



A/(p 2 ,g 2 ): 



p2 _ g2 



(19) 
(20) 

(21) 



where f = A,B. A feature of Eq. fll9|) is that the vertex is completely determined by 
the renormalised dressed-quark propagator. In Landau gauge the quantitative effect of 
modifications, such as that canvassed in Ref. [IS], is small and can be compensated for by 
small changes in the parameters that characterise a given model study f2(| . 

In the chiral limit (P 2 = 0) using Eqs. ([!]) and (|7|), the divergence of the AVV vertex is 



P p T mv (h, k 2 ) = H^ih, k 2 ) + /° T%,{k u k 2 ) 
where the direct contribution from the axial-vector vertex is 



RL(h,k 2 ) := —N c / tr DF 



T 



S{qi)ls-[l-PF R {q; 0) 



(22) 



(23) 



+7 ■ qq ■ PG R (q; 0) + q^P u H R (q; 0)J S(q 2 ) ^(g 2) qi2)S(q 12 ) iTl(q 12 , q x ) 
and that from the pion bound state is 



TUkx,k 2 ) := N c [ A tr DF 
J a 



S( gi ) l5 T 3 (iE„(q;0)- 1 -PF 7r (q;0) 
-J-qq- PG n (q; 0) - q^P v H n (q; 0)^ S(q 2 ) iE^(g 2 , q l2 )S(q 12 ) iTl{q 12 , q x ) 

Using Eqs. ©-(0), Eq. © simplifies: 

P P T Pflu (h, k 2 ) = RUkx, k 2 ) + f°T 3 (h, k 2 ) ■ 



(24) 



R 3 (h,k 2 ) :=-N c tr DF 



S(qi)l,-[l-PA (q 2 ) 



+ 1 -qq-P 2A' (q 2 ) S(q 2 ) iTVq 2} q 12 )S(q 12 ) ^(fe, qi] 



f 3 (h,k 2 ):=N c tr DF 



S(qi) J 5 r 3 iE n (q; 0) S(q 2 ) iTVq 2 , q 12 )S(q 12 ) iT2(q 12 , q x ) 



(25) 
(26) 

(27) 



5 



Now using Eqs. (P and flT9|) in Eq. (p6|) yields 



where (s := g 2 , A' := 4-A (s), etc 



o:. 



7T 



/•oo , 

ft(0) = jf rfs s 2 A 2 a v [A c 



a 



\2 



0/ o 



sA' a v + B' a c s 



0. 



(28) 



(29) 
(30) 



The last line follows because, using Eq. (|j) to eliminate a v and in favour of Aq and 
Bo, the integrand is identically zero. Hence the pseudovector components of the neutral- 
pion Bethe-Salpeter amplitude combine with the regular pieces of the axial-vector vertex to 
generate that part of the AVV vertex which is conserved. 

To reveal the anomalous contribution to the divergence, consider Eq. (27), in which using 
Eqs. (|) and (|IJ) yields 



71 



-fiupcr 



h P k 2a T{0) 



T(0) =j Q dss K A a° v (A [a° v a° s + s (< a° s - a° v a° s ') + s a° v (A' or° s -B' a 

Now, introducing C{s) := B (s) 2 /[sA (s) 2 ], Eq. (|2|) simplifies 

K(s;0) C'{s) 



T(P 2 = 0) = - 
which, using Eq. (|9|), yields 

ft T(P 2 = 0) = / 
Jo 

so that, in the chiral limit, 

Pp Tppu(ki, k 2 ) = 



ds s 



B (s) [1 + C(s)]3 



7T 



^ p,v p(jk\pk 2(J 



1 
2 



(31) 
(32) 

(33) 

(34) 
(35) 



Hence the pseudoscalar piece of the neutral-pion Bethe-Salpeter amplitude provides the 
only nonzero contribution to the divergence of the AVV amplitude. This contribution is 
just that identified with the "triangle anomaly", and the result is independent of detailed 
information about T n and S(p). It follows straightforwardly from Eqs. fl3T|) and (|3~4[) that 



3 2 
167T 7T 2 1 } 



64tt \7r/o 



(36) 



We emphasise that in obtaining the results in this section DCSB was crucial, since it 
originates and is manifest in a nonzero value of Bq, in the identity between Bq and E n , and 
in the other identities: Eqs. (PD[)-(P^). 

Our derivation is a generalisation of that in Ref. |14| and, to make it simple, particular 
care was taken in choosing the momentum routing in Eq. (|T7|). This was necessary because 
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it is impossible to simultaneously preserve the vector and axial vector Ward-Takahashi iden- 
tities for triangle diagrams in field theories with axial currents that are bilinear in fermion 
fields. This choice of variables ensures the preservation of the vector Ward-Takahashi iden- 
tity, which is tied to electromagnetic current conservation. With another choice of variables, 
surface terms arise that modify the value of 1Z(P 2 = 0). However, these are always elimi- 
nated by subtraction in any regularisation of the theory that ensures electromagnetic current 
conservation. flT5] 



III. ELECTROMAGNETIC PION FORM FACTOR 

As another example of the importance of IVs pseudovector components, we consider the 
electromagnetic pion form factor, calculated in the renormalised impulse approximation: 

(pi+P2)»F 7r (q 2 ) :=A> l5 p 2 ) (37) 

2N C r d^k r- i 
= -tfi J J^^ tT D[Gn(k]-p 2 )S(k ++ )iTl(k ++ ,k + ^)S(k + ^)g n (k-q/2]p 1 )S(k^^)\, 

k a /3 '■= k + api/2 + (3q/2 and p 2 := pi + q. Again, no renormalisation constants appear 
explicitly in Eq. ( |3"7| ) because the renormalised dressed-quark-photon vertex, T^, satisfies 
the vector Ward-Takahashi identity, Eq. ([18]). This also ensures current conservation: 

(p 1 -p 2 ) M A M (p 1 ,p 2 ) = 0. (38) 

We note that from the normalisation condition for Q n , Eq. ([13]), and Eqs. (0) and (]37 ) 



F(q 2 = 0) = 1 (39) 
if, and only if, one employs a truncation in which K is independent of P. One such scheme 



is the rainbow-ladder truncation of Ref. [11 . 



A. Quark propagator 

To calculate F^{q 2 ) we employ an algebraic parametrisation of the renormalised dressed- 
quark propagator that efficiently characterises many essential and robust elements of the 
solutions obtained in studies of the quark DSE. This defines Eq. (|3"T| ) directly Vpf,p|; in 
particular at the pion mass shell.0 We introduce the dimensionless functions: as{x) := 
\<Js(p 2 ), &v(%) '■= ^ 2cr v{p 2 ), where p 2 = X 2 x, X is a mass-scale, with 

<7 S (0 = 2m^(2(e + m 2 )) + T{h x OHh £) (&o + h 2 F{e 0) , (40) 
_ 2{i + m 2 )-l+e- 2 ^ 2 ) 

^ = WT^f — (41) 



The procedure actually employed in Ref. |16[ can, at best, only reproduce our results. 



7 



and J-{y) ■= [1— exp(— y)]/y. This five-parameter algebraic form, where m is the u/d current- 
quark mass, combines the effects of confinement^ and DCSB with free-particle behaviour at 
large, spacelike p 2 .f] 



The chiral limit vacuum quark condensate in QCD is |10| . |TT 

3 r M 2 

- (qq)l := Jirn^ Z 4 (/i 2 ,M 2 ) — J q ds s a° s (s) , (42) 

where at one-loop order Z 4 (/i 2 ,M 2 ) = [a(M 2 )/a(/i 2 )] 7m ^ 1+ ^ 3 ^, with £ the covariant-gauge 
fixing parameter (£ = specifies Landau gauge) and 7 m = 12/(33 — 2Nf) the gauge- 
independent anomalous mass dimension. The ^-dependence of Z 4 (/z 2 ,M 2 ) is just that re- 
quired to ensure that ° is gauge independent. The parametrisation of Eq. (|40| ) is a model 
that corresponds to the replacement 7 m — > 1 in Landau gauge, in which case Eq. fl42D yields 



This is the signature of DCSB in the model parametrisation and we calculate the pion mass 
from 

m 2 / 2 = 2m(gg)? GeV . (44) 
When all the components of T n are retained, Eq. (f44|) yields an approximation to the pion 



mass found in a solution of the Bethe-Salpeter equation that is accurate to 2% ]TT . 

The model parameters are fixed by requiring a good description of a range of pion 
observables. This procedure explores our hypothesis that the bulk of pion observables can 
be understood as the result of nonperturbative dressing of the quark and gluon propagators. 



B. Pion Bethe-Salpeter Amplitude 

The Chebyshev moments of the scalar functions in T n (k; P) are, for example, 

Ei(k 2 ;P 2 ) :=- I" dp sm 2 pU i (cosp)E 7T (k;P) (45) 

7T JO 



with k ■ P := Vk 2 P 2 cos/5, where Ui(z) is a Chebyshev polynomial of the second kind. At 
large- A; 2 , independent of assumptions about the form of K, one has fTI 



1 The representation of S(p) as an entire function is motivated by the algebraic solutions of Eq. (|3|) 



in Refs. [17]. The concomitant absence of a Lehmann representation is a sufficient condition for 
confinement. 
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5 At large-p 2 : <Jv{p 2 ) ~ 1/p 2 an d crs(p 2 ) ~ m/p 2 . The parametrisation therefore does not incor- 
porate the additional lnp 2 -suppression characteristic of QCD. It is a useful but not necessary sim- 
plification, which introduces model artefacts that are easily identified and accounted for. e = 0.01 
is introduced only to decouple the large- and intermediate-p 2 domains. 



S 



E° n (k 2 ;P 2 )^ -(qq)l 2 ^l. (46) 

F®(k 2 ;P 2 ), k 2 G°(/c 2 ; P 2 ) and k 2 H®(k 2 ; P 2 ) have precisely the same behaviour; i.e., the 
asymptotic momentum-dependence of all these functions is identical to that of B (k 2 ). 
This makes manifest the "hard-gluon" contribution to F n (q 2 ) in Eq. ( p7|) . Further, in 
an asymptotically free theory, where a well constructed rainbow-ladder truncation yields 



model-independent results at large- /r [ I I 



k 2 G»(k 2 ; P 2 ) = 2F»(k 2 ; P 2 ) , k 2 > M uv , (47) 

with M uv := 10 Aqcd- 

In a model exemplar used in Ref. [ll|] the zeroth Chebyshev moments provided results 
for and f n that were indistinguishable from those obtained with the full solution. Also 

~ and hence it was quantitatively unimportant in the calculation of m n and f^. We 
expect that these results are not specific to that particular model. In the latter case because 
the right-hand-side of Eq. (|1~2|) is zero and hence in general there is no "seed" term for H^. 



These observations, combined with Eqs. (p[)-([T2|), motivate a model for r„-: 

E 7r (k;P) = ^-B (k 2 ) (48) 

with F n (k;P) = E n (k; P)/ (110 U), G^k-P) = 2F n (k;p)/[k 2 + M 2 y ] and H v (k]P) = 0. 
The relative magnitude of these functions at large k 2 is chosen to reproduce the numerical 
results of Ref. [TT|. 



C. Results 



We determined the model parameters by optimising a least-squares fit to r n and 
(<?<?)iGeV' an d a selection of pion form factor data on the domain q 2 6 [0,4] GeV 2 . The 
procedure does not yield a unique parameter set with, for example, the two sets: 

A(GeV) m b b x b 2 b 3 
A 0.473 0.0127 0.329 1.51 0.429 0.430, (49) 
B 0.484 0.0125 0.314 1.63 0.445 0.405, 

providing equally good fits, as illustrated in Table ||] There is a domain of parameter sets 
that satisfy our fitting criterion and they are distinguished only by the calculated magnitude 



6 The quark propagator obtained with these parameter values is pointwise little different to that 
obtained in Ref. jlij] . One gauge of this is the value of the Euclidean constituent quark mass; i.e., 
the solution of p 2 = M(p 2 ) 2 . Here M% d = 0.32 GeV whereas M^ /d = 0.30 GeV in Ref. 0. It is 
also qualitatively similar to the numerical solution obtained in Ref. [jnj], where M^ d = 0.56 GeV. 
Indeed, our results are not sensitive to details of the fitting function: fitting with different confining, 
algebraic forms yields S(p) that is pointwise little changed, and the same results for observables. 
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of the pion form factor at large-g 2 . The two sets in Eq. (49) delimit reasonable boundaries 



and illustrate the model dependence in our result. With all parameter sets in the acceptable 
domain, Eq. ([TJ|) is satisfied exactly in the chiral limit, in which case we obtain /° = 
0.090 GeV, while at the fitted value of m, N w /f w = 0.97. 

In our calculation f n r n is 20% too small. This discrepancy cannot be reduced in impulse 
approximation because the nonanalytic contributions to the dressed-quark-photon vertex 
associated with n-ir rescattering and the tail of the p-meson resonance are ignored M. It can 
only be eliminated if these contributions are included. We have thus identified a constraint 
on realistic, impulse approximation calculations: they should not reproduce the experimental 
value of f n r w to better than ~ 20%, otherwise the model employed has unphysical degrees- 
of-freedom. 

Our calculated pion form factor is compared with available data in Figs. [1] and |2|. It 
is also compared with the result obtained in Ref. fll4" |, wherein the calculation is identical 



except thai the pseudovector components of the pion were neglected. Figure |I| shows a small, 
systematic discrepancy between both calculations and the data at low q 2 , which is due to 
the underestimate of rv in impulse approximation.^ The results obtained with or without 
the pseudovector components of the pion Bethe-Salpeter amplitude are quantitatively the 
same, which indicates that the pseudoscalar component, E n , is dominant in this domain. 

The increasing uncertainty in the experimental data at intermediate q 2 is apparent in 
Fig. H, as is the difference between the results calculated with or without the pseudovector 
components of the pion Bethe-Salpeter amplitude. These components provide the dominant 
contribution to F n (q 2 ) at large pion energy because of the multiplicative factors: 7 • P and 
7 • kk ■ P, which contribute an additional power of q 2 in the numerator of those terms 
involving F 2 , FG and G 2 relative to those proportional to E. Using the method of Ref. [14 
and the model-independent asymptotic behaviour indicated by Eq. (00) we find 



2 , . n(r) (-(qq){ 



7 4 /<?y ; (50) 

i.e., q 2 F n (q 2 ) « const., up to calculable In g 2 -corrections. If the pseudovector components 
of T n are neglected, the additional numerator factor of q 2 is missing and one obtains fl4[ 
q A F 1T (q 2 ) ~ const. 

In our model the behaviour identified in Eq. (|50|) becomes apparent at q 2 > 2 My V . This 
is the domain on which the asymptotic behaviour indicated by Eq. fl4"rJ| ) is manifest. Our 
calculated results, obtained with the two sets of parameters in Eq. (^), illustrate the model 
dependent uncertainty: 

q 2 FJq 2 ) ~ 0.12 -0.19 GeV 2 . (51) 



This uncertainty arises primarily because the model allows for a change in one parameter 

lA. 

'2 ! 



to be compensated by a change in another. In our example: 6f > b§ but b$ + 6f = 6q + ^2 ' 



7 Just as in the present calculation, f^r^ = 0.25 in Ref. [14]. However, the mass-scale is fixed so 



that f n 
larger. 



0.084, which is why this result appears to agree better with the data at small-g 



2. 



is 
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and &i 63 = bf 6f. This allows an equally good fit to low-energy properties but alters the 
intermediate-g 2 behaviour of F n (q 2 ). In a solution of Eq. (||) these coefficients of the 1/p 4 
and 1/p 6 terms are correlated and such compensations cannot occur. 

As a comparison, evaluating the leading-order perturbative-QCD result with the 
asymptotic quark distribution amplitude: <fias(x) '■= y/^fnx(l — x), yields q 2 F n (q 2 ) = 
I6irf 2 a(q 2 ) ~ 0.15GeV 2 , assuming a value of a(q 2 ~ lOGeV 2 ) ~ 0.3. However, the 
perturbative analysis neglects the anomalous dimension accompanying condensate forma- 
tion.0 



IV. CONCLUSIONS 



Using the Dyson- Schwinger equations it is straightforward to show that, as a conse- 
quence of the dynamical chiral symmetry breaking mechanism, the pion is a nearly-massless, 
pseudoscalar, quark- ant iquark bound state |]TU| , |rT| . As a corollary, the complete pion Bethe- 
Salpeter amplitude necessarily contains pseudovector and pseudotensor components, which 
are always qualitatively important. In model studies, the quantitative effect of these compo- 
nents can be obscured in the calculation of many pion observables; i.e., within a judiciously 
constructed framework, applied at low- to intermediate-energy, their effect can be absorbed 
into the values of the model parameters However, they are crucial to a proper realisa- 

tion of anomalous current divergences, crucial to obtaining a uniformly accurate connection 
between the low- and high-energy domains, and they provide the dominant contribution to 
the electromagnetic pion form factor at q 2 > 10 GeV 2 . 
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TABLES 



Calculated Experiment 



U 0.092 GeV 0.092 

(-(^)iGev) 1/3 0236 °- 236 ± °- 008 lH 



m u/d 0.006 0.008 ± 0.004 |22 

0.1387 0.1385 



r n 0.55 fm 0.663 ± 0.006 (2J 

tv/ w 0.25 (dimensionless) 0.310 ± 0.003 

TABLE I. A comparison between our calculated values of low-energy pion observables and 
experiment or, in the case of (— (q^igcv) 1 ^ 3 an d rn u / d , the values estimated using other theoretical 
tools. Each of the parameter sets in Eq. (|4^) yields the same calculated values. For consistency 
with Ref. [11 1, we use Aqcd = 0.234 GeV throughout. 
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FIGURES 
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FIG. 1. Calculated pion form factor compared with data at small q 2 . The data are from 

(circles). The solid line is the result obtained when the pseudovector 



Refs. 23] (crosses) and |2 
components of the pion Bethe-Salpeter amplitude are included, the dashed-line when they are 



neglected [14]. On the scale of this figure, both parameter sets in Eq. ([49]) yield the same calculated 
result. 
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FIG. 2. Calculated pion form factor compared with the largest q 2 data available: diamonds 
and circles 



Ref. and circles - Ref. [25]. The solid lines are the results obtained when the pseudovector 
components of the pion Bethe-Salpeter amplitude are included (lower line - set A in Eq. p^); upper 
line - set B), the dashed-line when they are neg lected 0. 
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